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ABSTRACT 

We present results from the first solar full-disk F 10.7 (the radio flux at 10.7 cm, 2.8 GHz) 
image taken with the S-band receivers on the recently upgraded Karl G. Jansky Very Large 
Array (VLA) in order to assess the relationship between the F 10.7 index and solar extreme ultra¬ 
violet (EUV) emission. To identify the sources of the observed 2.8 GHz emission, we calculate 
differential emission measures (DEMs) from EUV images collected by the Atmospheric Imaging 
Assembly (AIA) and use them to predict the bremsstrahlung component of the radio emission. 
By comparing the bremsstrahlung prediction and radio observation we find that 8.1 ±0.5% of the 
variable component of the F 10.7 flux is associated with the gyroresonance emission mechanism. 
Additionally, we identify optical depth effects on the radio limb which may complicate the use 
of F 10.7 time series as an EUV proxy. Our analysis is consistent with a coronal iron abundance 
that is 4 times the photospheric level. 

Subject headings: Sun: abundances, Sun: corona, Sun: radio radiation, Sun: solar-terrestrial relations, 
Sun: UV radiation 


1. Introduction 

Microwave emission from the sun, specifically 
the F 10.7 index (the solar radio flux at 10.7 cm, 2.8 
GHz), has long been known to correlate with so- 


scales of days and longer (Covington 1948 

, 1951: 

Covington et al. 19551 

Vats et al. 

1998: 

Foukal 

1998; Tapping et al. 200311 and has been measured 

daily since 1947 (Covington 19691; Tapping 1987). 


This correlation and the transparency of th e at 


mosphere to microwave signals (Tapping 120131) 


has led to the use of F 10.7 as a proxy measure¬ 
ment for solar EUV irradiance, which heats and 
ionizes the Earth’s atmosphere but cannot be ob¬ 
served from the ground. The observed F 10.7 signal 
was one of the inputs for the original terrestrial 
i onos phe ric and ther mospheric modelling efforts 
( Bhatnagar fc Mitra 19661 Jacchial 1~1 97 ill and re¬ 


mains one of the primary model inputs today, 


even wi th the availab i lity o f direct EUV obser¬ 
vations (ITobiska et al. 200811 . F 10.7 is often pre¬ 
ferred over other proxies of solar activity such as 
Sunspot Number (SSN), the Mgll core-to-wing 
index, Lya irradiance, etc. due to its high de¬ 
gree of correlation with solar EUV output which 
results from the fact that the main variable com¬ 
ponents of both originate in the same coronal 
plasma (ISwarup et al.lll963lf . 

As early as Kundu (.1959) it was understood 
that thermal microwave sources could be split 
into three classes: a low intensity backgroun d 
originating from the quiet Sun ( Martvn 19481) , 
a moderate intensity sign al seen in and a round 
plage and active regions ( Cov ington 119471 ). and 
a high intensity componen t co mmonly associated 
with active region cores ( Piddin gton fe Minnett 
Il95ll) . The background component (here taken 
to be 65.2 ± 2.0 sfu of the F 10.7 flux) can be ex¬ 
plained with a uniform optically thick chrorno- 
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gated the variability of long duration coronal mi¬ 
crowave (primarily Fin?) timeseries for use as 


proxies of EUV (Covington 1951, 19691; Tapping 

1987; Wilson et al. 1987f Lean & Brueckneij 

1989 


Bouwer 1992; 

Foukal 19981; Parker et al. 

1998 


Vats et al. 1998 

Tobiska et al. 

2008: Dudok de Wit et al. 

2009: Maruvam 

a. 2010: Svalgaard & Hudson 

2010 


Chen et al. 201 

Johnson 2011; Maruvama 

2011 


Hennev et al. 20121; Denget al. 

20131) and total so 


lar irradiance (Tapping et al. 2003 Frohlich 2009t 

Tapping & Valdes 2011) but 

there is still 

con 



siderable debate as to the exact source of the 
time-variable component. Some studies have 
argued that the variable microwave component 
is optically thin bremsstrahlung emissio n origi¬ 
nating from the plae 
Tapping fc DeTrace'v 


age regi ons (Fclli et al. 198 li 
•vi ll990t Tapping et al.l 12003 ) 


while others suggested that it is primarily gyrores- 
onance emission from the str ong mag netic fields 
in active region cores ( Schmahl fo Kund~ul Il995l 


1998ll . As discussed below (in Q, bremsstrahlung 


emission is closely related to EUV emission, while 
gyroresonance emissi on is_not. A recent study by 


Dudok de Wit et alJ (20141) suggests that the gy¬ 


roresonance component can “account for 90% of 
the rotational variability in the F 10.7 index.” 

This ambiguity in the source of the signal vari¬ 
ability is best resolved through imaging when the 
individual sources are resolved. Many such studies 
have been performed throughout the m icrowave 


1963 


1991 



There has been less work specifically at 2.8 
GHz, jmd the best i magin g to date was that 
of Saint-Hilaire et al. ( 2012h who used the Allen 
Telescope Array to observe the full Sun between 
1.43 and 6 GHz (including F 10 . 7 ). They used the 
emission spectra and polarization to identify gy¬ 
roresonance sources with a spatial resolution at 2.8 
GHz of about V. However, due to a lack of avail¬ 
able instrumentation, imaging of F 10.7 with reso¬ 
lution better than one arcminute did not become 
possible until the upgrade to the Karl G. Jansky 


Very Large Array (VLA). ISelliorst et al. ( 2014 ) 
recently completed a study of the statistical prop¬ 
erties of spatially resolved active regions observed 
at 17 GHz in order to identify bremsstrahlung 
and gyroresonance emission and noted that a sim¬ 
ilar analysis is needed to fully understand F 10 . 7 . 
While high-resolution radio studies are necessary 
to identify gyroresonance regions and their contri¬ 
bution to F 10 . 7 , they are not in general sufficient 
to determ ine the magnitude of t he gyroresonance 
emission ( White &; Kundu 19971) and some inde¬ 
pendent estimate of one of the components of the 
radio signal is required. 

The bremsstrahlung component of the flux can 
be extrapolated without direct radio observation if 
the differential emission measure (DEM, i.e., the 
distribution of plasma density with temperature 
in the corona) is known. DEMs can be calcu¬ 
lated from sets of optically thin observations, in¬ 
cluding EUV imaging and spectroscopy, as long 
as they are sensitive to a range of coronal tem¬ 
peratures wide enough t o sa mple th e do minant 
coronal plasma. Landi fc Chiuderi Drago ( 2003 . 
20081) constructed solar-minimum DEMs using 


both EUV and radio observations to constrain the 
plasma structure from the chromosphere, which 
has significantly fewer EUV emission lines, all the 
way through the corona. With current instrumen¬ 
tation, coronal DEMs are most commonly con¬ 
structed from EUV data because of the high spa¬ 
tial, temporal, and spectral resolution of obser¬ 
vations, especially in the jSolar Dynam ics Obser¬ 
vatory (SDO) era (P esnel l et al. 120111 ). Compar¬ 
isons between radio observations and the predicted 
bremsstrahlung component based on DEMs com¬ 
puted from EUV images have b een made previ¬ 


ously, e.g., bv lWhite et al. ( 2000h and Zhang et al 
(2001). EUV emission is dominated by specific 


atomic emission lines while the bulk thermal emis¬ 
sion in the radio is generated by free electrons and 
therefore represents the density of fully-ionized 
hydrogen in the corona. This means that the 
elemental abundance is a (presumably) constant 
scaling factor relating the element-specific DEM 
and the radio measurement. If the abundance is 
known, then the DEM can be used to predict the 
optically thin bremsstrahlung emission fro m the 
plasma observed in the EUV. IWhite et al. ( 2000h 
used this technique to measure the Fe abundance 
in an active region. 
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In this paper we analyze the first full-disk im¬ 
age of F 10.7 emission acquired with the VLA, 
which is the highest spatial resolution 2.8 GHz 
image to date, and compare it with spatially- 
resolved EUV images. The prediction of the 
bremsstrahlung emission along with the polariza¬ 
tion signal in the radio are used to identify gy- 
roresonance sources in the radio image and de¬ 
termine the total gyroresonance component of the 
F 10.7 flux. In 0 the physical connections be¬ 
tween the various EUV and radio emission mech¬ 
anisms are discussed, with details of the observa¬ 
tions given in fJ5J The calculation of the DEM is 
described in 0 and the subsequent comparison of 
the bremsstrahlung prediction and the radio ob¬ 
servation is explained in (J5] |J6] includes a discus¬ 
sion of the results and implications, and we con¬ 
clude and suggest future additions to this work in 

El 


2. Physical Background 


2.1. Radio Emission 

There are only two identified mechanisms that 
produce microwave radio emission from the non- 
flaring sun: bremsstrahlung and gyroresonance 
( Kundul 19651) . 

Thermal bremsstrahlung results from th e co lli- 


sional interaction of electrons and ions (IWild et al 
1963). For optically thin bremsstrahlung emis¬ 
sion in coronal conditions the radio flux density 
(in units of erg cm -2 s -1 Hz -1 ) is given by: 


f v = 9.78 x 10 


3 2fcB 


1+4 


Nile 

iV H 


T-°' 5 DEM(T)G(T) dT dfl (1) 


where = 1.38 x 10 16 g cm 2 s 2 K 1 is Boltz¬ 
mann’s constant, c = 3 x IQ 1 0 cm s -1 is the speed 


of light, Ane/An = 0.085 (Asplund et al. l2009h 
is the number (or number density) ratio of He¬ 
lium to Hydrogen in the emitting medium, T is the 
temperature in Kelvin, G(T) = 24.5 + In (T/v) is 
the Gaunt factor where v is the freq uency i n Hz , 
clD is the solid angle of the source ( Dulkl Il985l) 
and DEM(T) is the integral along the line of sight 
through the corona of d(n e nH)/dT. The depen¬ 
dence of the flux on n 2 T -05 means that the op¬ 
tically thin bremsstrahlung flux is actually rela¬ 
tively insensitive to the temperature distribution 


and is much more sensitive to the plasma density 
(see ®. 

Thermal bremsstrahlung emission at microwave 
frequencies (1 — 30 GHz) generally becomes opti¬ 
cally thick in the chromosphere because of both 
the increased density and decreased temperature. 
The altitude (and therefore temperature and den¬ 
sity) at which this optically thick boundary oc¬ 
curs is a strong function of frequency v since 
bremsstrahlung opacity varies as v~ 2 n 2 T~ 0 - 5 , 
with higher frequencies penetrating deeper into 
the chromosphere. This leads to a frequency de¬ 
pendence in the observed height of the solar limb 
at microwave frequencies, with the apparent size 
of the s olar disk decreas ing with increasing fre¬ 
quency ( Furst et~ahlll979l) . In active regions with 
very high density this optically thick boundary 
can also occur in the corona at low microwave fre¬ 
quencies, increasing the observed brightness tem¬ 
perature dramatically and blocking observation of 
the lower atmosphere. It is common for active 
regions to be optically thick in the corona due 
to bremsst rahlung at 1 .4 GHz but optically thin 
at 5 GHz ( White Il999i) . putting F 10.7 (2.8 GHz) 
at an interesting transition frequency. Coronal 
bremsstrahlung emission is generally not strongly 
polarized, but magnetic fields do break the degen¬ 
eracy of collisional in teractions and produce weak 
circular polarization ( White fc Kundu 19971 ). 

Gyroresonance emission arises from the acceler¬ 
ation of electrons as they spiral around magnetic 
field lines. At coronal temperatures even ther¬ 
mal electrons have weakly relativistic velocities 
and produce opacity not just at the gyrofrequency , 
but also at low order harmonics ( Wild et all 19631 ). 
Depending primarily on the magnetic field orien¬ 
tation and the polarization mode, this emission 
becomes optically thick in the s = 1, 2, 3, or 4 
harmonic of the gyrofrequency: 


v B = 2.80 B [MHz] 


( 2 ) 


where B is the magnitude o f the magnetic field 
in G ( White fc Kundul Il997 ). This means F 10.7 
gyroresonance observations at 2.8 GHz come 
from thin surfaces with constant magnetic field 
strengths of B = 10 3 /s = 1000, 500, 333, and 250 
G for harmonics s = 1, 2, 3, and 4 respectively. 

The characteristic motion associated with the 
gyroresonance process naturally causes the emis¬ 
sion to be highly circularly polarized (because any 
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intrinsic linear polarization is wiped out by Fara¬ 
day rotation in the solar atmosphere). Emission 
with polarization in the sense of an electron spi¬ 
ralling around the field is called the extraordinary 
or ir-mode, while polarization with the opposite 
sense of rotation is called the ordinary or o-mode. 
Electrons couple much more strongly to the x- 
mode than the o-mode because of the shared sense 
of rotation. Consequently, the ir-mode generally 
has larger opacity and becomes optically thick in 
higher (harmonic and altitude) gyroresonant lay¬ 
ers. The generally positive temperature gradient 
in the lower corona means that the ir-mode then 
has a higher brightness temperature causing an 
observed net circular polarization from gyroreso- 
nance sources in the sense of the local ir-mode. 


2.2. EUV Emission 


While there is high energy X-ray bremsstra hlung 
continuum emission ( Craig fc Brown] 1197611 . the 
majority of the energy output from the non¬ 
flaring corona comes from collisionally excited 
atomic emission lin es, predo minantly observed in 
the EUV ( Golub fc Pasachoflboiotl . At coronal 
temperatures the EUV spectrum is dominated by 
emission lines with strengths dependant on a num¬ 
ber of characteristics of the bulk plasma and the 
individual emitting atom. These include but are 
not limited to: the local plasma density (which 
influences the collision rate), the local electron 
temperature (which determines the energy of the 
collisions), the ionization and excitation state of 
the atom (which restrict the available transitions), 
and the oscillator strengths of the available tran¬ 
sitions (which determine both the probabilities of 
each emission as well as their respective energy 
spectra). However, assuming the atomic details 
are known, the EUV emission properties of a bulk 
plasma in ionization equilibrium are completely 
determined by the relative elemental abundances 
and t he electron density an d temperature distribu¬ 
tion ( Craig fc Brown 19761 ). The temperature de¬ 
pendence of individual emission lines means that 
observing specific regions of the EUV spectrum 
highlights very different coronal features. One 
commonality among these coronal EUV obser¬ 
vations is that they remain optically thin in all 
non-flaring conditions, and therefore all coronal 
plasma above the much denser and optically thick 
chromosphere is visible. 


2.3. DEM Connection 

As has been noted, both radio bremsstrahlung 
continuum and EUV emission line strengths are 
at least partially dependent on the plasma den¬ 
sity and temperature distribution which can be de¬ 
scribed by the DEM, and therefore the two can be 
compared. We compute the DEM (J3} using the 
higher spatial, temporal, and spectral resolution of 
EUV observations and then use it to predict the 
optically thin bremsstrahlung radio emission from 
the same plasma. To first order, the difference be¬ 
tween the optically thin bremsstrahlung prediction 
and a radio measurement should reflect the opti¬ 
cal depth effects in the bremsstrahlung emission 
and/or the presence of gyroresonance components. 
Note that the detailed temperature structure of 
the DEM is less important for this application be¬ 
cause the comparison with the radio only depends 
on the (temperature-weighted) integral over the 
DEM as shown in equation Q] 

This comparison lacks the information needed 
to account for the altitude of the optically thick 
layer in the radio, which adds a complicating fac¬ 
tor. Both bremsstrahlung emission from dense ac¬ 
tive regions and gyroresonance emission can gener¬ 
ate optically thick layers in the corona that block 
radio emission from any plasma below that layer. 
The EUV lines used to compute the DEM remain 
optically thin all the way to the chromosphere, 
therefore the EUV will generally observe more 
plasma than the radio and may lead to a rela¬ 
tive overestimate of the observable optically thin 
bremsstrahlung emission. This effect is difficult to 
assess and we discuss it further in H6.4I 


3. Data 


The observations for this analysis were taken 
during the rising phase of solar cycle 24 on 2011 
December 9, between 15 and 23 UT. During this 
period there was moderate solar activity on the 
earthward hemisphere and the Fm .7 index was 


143.5 ±1.2 sfu ( iTapping fc Charroisl 1994 1. but 
there were no recorded solar storms of any kind. 
The observed variability of coronal features was 
insignificant and occurred mostly on scales below 
the resolution of the radio observations. 
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3.1. Karl G. Jansky Very Large Array 

At the time of these observations, the Karl 
G. Jansky Very Large Array (VLA, operated by 
the National Radio Astronomy Observatory) was 
in the process of being upgraded, and data were 
taken as a shared-risk project. As such, only 17 
of the VLA’s nominal 27 available antennae were 
equipped with the “S-band” feeds (2 — 4 GHz) 
needed to observe F 10 . 7 . The VLA was in its most 
compact (“D”) array configuration, appropriate 
for recovering the flux of large-scale sources in the 
solar atmosphere. The reduced number of anten¬ 
nae decreased the resulting image quality signifi¬ 
cantly, because of both the reduced coverage of the 
u-v plane, as well as the 37% decrease in collect¬ 
ing area. For this analysis, eight 2 MHz channels 
centred on 2.783 GHz were summed for a total 
bandwidth of 16 MHz. Solar observations in S- 
band are taken with the additional nominally -20 
dB attenuators in the signal path and these add 
phase changes that are corrected using indepen¬ 
dently measured values of the delays (Bin Chen, 
private communication). Unfortunately, the mea¬ 
surements of the primary flux calibration source 
were corrupted, preventing independent measure¬ 
ment of the amplitude changes due to the atten¬ 
uators. Consequently, the solar fluxes were cali¬ 
brated assuming exactly 20 dB of attenuation and 
a nominal flux for the secondary calibrator. We 
estimate that the VLA solar fluxes therefore have 
an uncertainty of order 20 %. 

The field of view of a single VLA pointing with 
the S-band receivers is nominally 15' and there¬ 
fore mosaicking is required to image the entire 30' 
diameter solar disk. A honeycomb pattern mosaic 
was used, with a single pointing at disk center sur¬ 
rounded by six fields designed to overlap by half a 
beam width. The center of each field was tracked 
over the course of the eight hour observation, tak¬ 
ing into account solar differential rotation. This 
led to some slight feature smearing near the edges 
of each field, but the magnitude was well below the 
final 25" resolution of the observation and there 
appeared to be no effect on the final mosaic im¬ 
age. Each field was calibrated separately and then 
the fields were imaged jointly in a single map us¬ 
ing maximum entropy deconvolution. The images 
were restored with a spatial resolution of 25". Due 
to the limited field of view, as well as the restric¬ 
tion caused by the minimum baseline, the obser¬ 


vations were insensitive to emission on the scale 
of the solar disk. Attempts to restore this compo¬ 
nent using a default disk of the right dimension in 
the deconvolution process were unsuccessful (be¬ 
cause they failed to produce a mostly flat disk as 
seen in the Allen T elesc ope Array observation of 
Saint-Hilaire et al. 2012 :). and we therefore do not 
address the spatial distribution of the large-scale 
emission here. This large-scale component will be 
analyzed with a subsequent data set acquired us¬ 
ing more antennae and a larger mosaic pattern. 


3.2. Nobeyama Radioheliograph 

The 17 GHz data fro m the Nobeyama Ra ¬ 


dioheliograph (NoRH; Naka iima et al. 119941) 


were used to assess optical depth effects in the 
bremsstrahlung radio emission (the optical depth 
of bremsstrahlung at 17 GHz is 37 times smaller 
than at 2.8 GHz) and as a calibration check. 
NoRH makes full disk images of the Sun at 17 
and 34 GHz every day between 23:00 and 06:30 
UT. Located in Japan, the dedicated solar array is 
unable to observe simultaneously with the VLA, 
and therefore the Nobeyama images collected just 
after the completion of the VLA observation were 
used. The data were mapped and calibrated using 
standard procedures: amplitude calibration as¬ 
sumed that the background disk component that 
generally dominates the total flux had a brightness 
temperature of 10 4 K, which is known to be con¬ 
sistent with well-calibrated flux monitoring at this 
frequency by th e Nob eyama polarimeters (NoRP; 
Nakaiima et allll9851 ). The magnetic field depen¬ 
dence of equation [2] means that 17 GHz is only 
sensitive to gyroresonance emission from strong 
magnetic fields, requiring a coronal field greater 
than 2000 G to observe the third harmonic. Coro¬ 
nal magnetic field strengths this high are typi¬ 
cally only seen in the case of very large active 
regions. Given the absence of such regions during 
the observation, the Nobeyama observations are 
expected to detect purely bremsstrahlung emis¬ 
sion from the solar atmosphere that is excess to 
the (chromospheric) background brightness tem¬ 
perature level of 10 4 K. The 17 GHz image was 
made by synthesizing and deconvolving data taken 
at 45 s intervals over a period of several hours, and 
rotating the final image back in time to match the 
VLA image. The spatial resolution in the final 17 
GHz image is 12", and the flat background disk 
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of 10 4 K (with a radius 1.0125 times the photo- 
spheric radius that fits the 17 GHz visibilities) is 
subtracted for the region analysis. 


3.3. Atmospheric Imaging Assembly 


The full disk EUV images used for this anal¬ 
ysis came f rom the Atmospher ic Imaging Assem¬ 
bly (AIA) dLemen et al.l l2012h aboard the Solar 
Dynamics Observatory satellite. All six coronal 
EUV channels (94A, 13lA, 17lA, 193A, 21lA, and 
335A) were used at one minute cadence over the 
course of the VLA observat ion. The point spread 
function corrections from Poduval et al. ( 2013h 
were applied to the level 1.5 images which were 
then summed after rotation to a common time at 
the midpoint of the VLA observation to produce 
longer integrations and increase the signal to noise. 
The blurring created by small scale feature fluctu¬ 
ations during this long integration had no effect on 
our results because the bremsstrahlung prediction 
resulting from the EUV data was convolved with 
a 25" Gaussian beam before analysis to match the 
resolution of the radio observations. The EUV im¬ 
age sequence was also used to check for any time 
variability that might affect our results (Figurc[l]). 
No major time variability is present and therefore 
the time-integrated EUV data is appropriate for 
comparison with the 8-hour VLA data set. 


4. The Differential Emission Measure 

The emission measure (EM) of hydrogen in the 
corona is defined as: 


EM = / n e (s)nn(s) ds 


(3) 


where n e is the electron number density, nn 
is the hydrogen number density, and l is the 
path len gth through the optic ally thin medium 
(Greenstein & Minkowski 19531). Since the corona 
is not isothermal, it is conventional to use the 
differential emission measure (DEM) which repre¬ 
sents the column integrated plasma density as a 
function of temperature: 


EM = J DEM(T) dT (4) 

with 

DEM(r} = / d "g"( (a) * (5) 

The observed intensity I\ at wavelength A is then 
given by: 


h = J R\(T)DEM(T) d T (6) 


3.4. Helioseismic and Magnetic Imager 

Photospheric line-of-sight magnetic field mea¬ 
surements were obtained by the Helioseismic and 
Magnetic Imager (HMI) (Scherrer et alj 120 11 j on 
the Solar Dynamics Observatory. These 4096 x 
4096 low-noise full sun magnetic field map s are 
produced every 5 minutes (|Schou et ah 2011 ). but 
only the single observation closest to 19:00 UT 
(the central time of the VLA observation) was 
used because, like the corona, the photosphere 
showed very little variability during the obser¬ 
vation window. While radio and EUV obser¬ 
vations are sensitive to the corona and chromo¬ 
sphere which lie megameters above the photo¬ 
sphere, strong photospheric magnetic fields indi¬ 
cate large active regions which extend into the 
corona where gyroresonance emission should be 
the strongest. HMI magnetograms are used to 
qualitatively connect radio polarization measure¬ 
ments to the coronal magnetic field. 


where R\(T) is the temperature response function 
of the instrument. This response function is de¬ 
pendent on both the technical details of the in¬ 
strument (wavelength resolution, filter passband 
shape, detector response, etc.) and the atomic 
physics of the emitting plasma (composition, tran¬ 
sition probabilities, occupation states, tempera¬ 
ture sensitivity, etc.). By observing multiple emis¬ 
sion lines with different temperature sensitivity 
(either as spectral lines or through narrow band 
imaging) it is possible to invert the system of 
equations (equation [6]) to determine the DEM of 
the source plasma. However, with the addition of 
measurement errors such as Poisson noise and any 
instrumental effects, this system becomes under¬ 
constrained, and a precise, analytic inversion is 
impossible. Additionally, even if a self-consistent 
solution can be found, its relation to the actual 
emitting plasma is dependent on the atomic pa¬ 
rameters in R\(T) which, in the EUV, may have 
errors of up to ~ 50% (jZanna et ah 201 1). The 
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Fig. 1. : Full disk solar images of left: the eight hour integrated 21lA AIA image and right: the standard 

deviation of each pixel in the 21lA image time series plotted as a percentage of the observed flux. Notice 
that the on-disk variation is small, especially in active regions which have high signal to noise. 


AIA response f unctions calculat ed using the CHI¬ 
ANTI package ( Dere et aJ1l20Q9h were used assum¬ 
ing coronal abundances 1 116.311 and the CHIAN TI 
default ionization balance (Boerner et al. 2012). 


4.1. Solution Technique 

This w ork uses the i n versio n method pre¬ 
sented in Plowman et al. 12 013 1. The code is 


based on the original procedure presented in 
Hannah fe Kontai] ( 20121) and uses a multi-step 
process to invert a set of coronal flux measure¬ 
ments and derive the best source DEM distribu¬ 
tion. The first pass involves a direct inversion of 
the input data using the instrument response func¬ 
tions themselves as the set of linearly- independent 
basis functions. This method naturally applies 
a minimum squared EM condition that prefers 
smooth solutions which tend to minimize unphys¬ 
ical negative EM contributions. A regularized so¬ 
lution with narrower basis functions is computed 
in order to remove the negative EM components 
entirely. This regularization is compared with 
the original solution and iteratively modified to 
slowly reduce the negative emission while ensur¬ 
ing that each step deviates by no more than the 
accepted y 2 threshold. Three different regular¬ 
ization strengths are applied, and only if they 
all fail to converge are any negative EM com¬ 


ponents allowed to remain in the final solution. 
The primary advantage of this solution method 
is its speed, computing a full resolution AIA 
(4096x4096 pixels) DEM image in abo ut one hour 


on a single processor workstation (jPlowman et al 


12013 ). The AIA images used here are dominated 
by the lines of Fe VIII, IX, XII, XIV, XVIII, and 
XXI, which together cover the temperature range 
log(T) = 5.6 — 7.0 corresponding to the bulk of 
coronal plasma. 

4.2. Full Disk DEM 

Full disk representations of the calculated 
DEMs are shown in Figure [2] The left image 
shows the total emission measure which was ob¬ 
tained by integrating the DEM in each pixel over 
the temperature axis. This EM dominates the 
contribution to the final bremsstrahlung predic¬ 
tion because it is linear in the integral in equation 
[T] and has over two orders of magnitude variation 
on the solar disk. On the right is an image of the 
emission-measure-weighted median temperature 
calculated from the derived full disk DEM. This 
map shows some small discontinuities along the 
east and north limb (which appear saturated in 
the image) but these regions play no role in the 
analysis as described in )}5] The temperatu re im- 
pacts the bremsstrahlung emission as T -0 - 5 (|Dulk 
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Il985h in the integral over temperature (Equation 
HD: and, because the median temperature varies by 
only half an order of magnitude across the solar 
disk, the temperature variation plays a relatively 
minor role in the predicted radio fluxes compared 
to the emission measure. 

5. Analysis 

An image of the expected coronal component 
of the bremsstrahlung emission was calculated us¬ 
ing the DEMs computed from the AIA images and 
equation [1] assuming a coronal abundance for Fe 
discussed in > 16.31 This image was computed at 
the full AIA resolution and then convolved with a 
25" full width half maximum Gaussian and down- 
binned to match the resolution of the radio ob¬ 
servations. The brightness temperature of the 
bremsstrahlung prediction and radio observation 
are shown in Figure [3] on the same scale, with 
individual regions outlined and labelled. Notice 
that while the overall morphologies agree quite 
well, the observation has large, high temperature 
emission peaks from the centres of many of the 
regions which are absent from the bremsstrahlung 
prediction. Additionally, the bremsstrahlung pre¬ 
diction has significantly more faint emission sur¬ 
rounding the active regions than was observed. It 
is important to emphasize that the prediction im¬ 
age is based on the AIA observations and there¬ 
fore it will vary from the observed optically thin 
bremsstrahlung wherever the EUV is sensitive to 
different plasma than the radio. 

The average F 10.7 flux measured by the official 
Solar Monitoring Program in Penticton, Canada, 
at 18, 20, and 22 UT on 2011 December 9 was 
143.5 ± 1.2 sfu. Subtracting the quiet sun back¬ 
ground of 67.2±2.1 sfu (the 65.2±2.0 sfu constant 
solar minimum level scaled to a Sun-Earth sepa¬ 
ration of 0.985 AU) leaves an observed variable 
F 10.7 component of 76.3 ± 2.4 sfu. The total F 10.7 
flux from the predicted bremsstrahlung image was 
77.7 ± 0.1 sfu, which compares well with the ob¬ 
served variable flux. This suggests that gyroreso- 
nance emission is not distorting the F 10.7 flux sig¬ 
nificantly since the optically thin bremsstrahlung 
component can account for all of the variable F 10.7 
on this day. However, this conclusion ignores sev¬ 
eral complicating details which are discussed in (J6] 


5.1. Region comparison 

The full disk images in Figure [3] cannot be com¬ 
pared quantitatively because of the failure to re¬ 
store the flat background disk (which is in any 
case absent from the EUV images) to the F 10.7 
radio image. Instead, individual regions on the 
Sun for which the imaging is reliable were ana¬ 
lyzed. The bremsstrahlung prediction, radio in¬ 
tensity, circular polarization, and photospheric 
line-of-sight magnetic field strength in each re¬ 
gion were compared in order to determine if gy- 
roresonance emission was present. These compar¬ 
isons for the regions with the largest observed po¬ 
larizations are shown in Figure 01 In region 3, 
the general morphologies of the observed and pre¬ 
dicted active regions agree very well despite the 
much higher observed radio brightness tempera¬ 
ture than predicted by the EUV data. Addition¬ 
ally, the circular polarization signal is large above 
the strong photospheric magnetic fields, as ex¬ 
pected. The remaining three regions shown in the 
figure each display varying levels of morphologi¬ 
cal deviation between the bremsstrahlung predic¬ 
tion and radio observation, suggesting that there 
are gyroresonance sources offset from the peak 
bremsstrahlung emission (because optical depth 
effects in the bremsstrahlung sources will not pro¬ 
duce offsets). In region 4, the peak observed 
emission is shifted to the southwest by about 
20 ", aligning with the observed polarization sig¬ 
nal. The observation of region 10 extends east and 
west farther than the predicted emission, align¬ 
ing with the polarization signal directly above the 
east-west oriented photospheric magnetic fields. 
Region 15 shows a strongly polarized radio source 
above an isolated sunspot with no corresponding 
predicted bremsstrahlung source. These regions 
show that the VLA resolution is sufficient to ex¬ 
tract active region details in both the intensity and 
the individual polarization channels, allowing for 
the identification of large gyroresonance sources 
simply from the polarization and morphological 
inconsistencies. 

In order to isolate the active region fluxes from 
any larger-scale background, identical background- 
subtraction approaches were used in both the 
bremsstrahlung prediction and the radio obser¬ 
vation to allow quantitative comparison between 
the data sets. This involved using the solar disk 
around the border of each region to estimate the 





Fig. 2. : Full disk solar images of left: the total emission measure and right: the emission-measure-weighted 

median temperature as calculated from the AIA images. 
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Fig. 3.—: Full disk solar images on the same brightness temperature scale of left: the optically thin 
bremsstrahlung prediction and right: the 2.783 GHz observation. Individually analyzed regions are boxed 
and numbered. 


disk emission within the region itself. This was 
done for concentric borders up to 3' outside each 
region, using the average result to calculate the re- 
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gion flux and the variation in the total region flux 
as a measure of the uncertainty. Note that while 
these uncertainties are quoted for the remainder of 
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Fig. 4.—: Individual regions of interest which each have a peak polarization brightness temperature greater 
than 10 5 K. The left pane shows the bremsstrahlung prediction as an inverted heat map with the radio 
intensity over-plotted with white contours at (8,16,32, 64,128, 256, 512,1024) x 10 3 K. The right pane shows 
the photospheric line-of-sight magnetic field in gray scale with the radio polarization brightness plotted with 
contours (blue for left hand polarization and red for right hand polarization) at ±(8,16, 32, 64,128, 256, 512) x 
10 3 K. The small black circle in the top right corner of the left pane is the approximate beam size and the 
white plus signs mark the location of peak observed radio intensity. The left pane also lists the peak 
brightness temperature of the radio observation (white) and the bremsstrahlung prediction (black) while the 
right pane lists the minimum (left hand) and maximum (right hand) polarization brightness temperatures. 


the paper, the systematic uncertainties associated 
with the DEM calculation and the uncertainty in 
the VLA flux calibration may be much larger (as 
discussed in ij6.ll) . The flux from each region is 
plotted in Figure [3] where the regions have been 
classified based on their maximum polarization 
brightness temperature and their proximity to the 
solar limb. 

We emphasize peak circular polarization 
brightness temperature rather than degree 
of polarization because the effective noise 
level in the latter is very high. This is due 
to both the noise level in Stokes I and V 
and our inability to fully restore the disk 


emission in total intensity which greatly af¬ 
fects the low intensity regions. In addi¬ 
tion, the large beam size at 2.8 GHz re¬ 
sults in smearing between any smaller gy- 
roresonance sources and the more extended 
bremsstrahlung emission. In some cases 
(notably region 10), oppositely polarized 
gyroresonance sources overlap within the 
25" beam size and cancel, resulting in ar¬ 
tificially low polarization. All the regions 
with strong circular polarization show de¬ 
grees of polarization greater than 30%, but 
this is likely an underestimate. 
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Fig. 5.—: Total flux observed in each region at left: 2.783 GHz and right: 17 GHz is plotted against the 
total predicted optically thin bremsstrahlung emission. The dotted lines indicate where the bremsstrahlung 
prediction equals the observed flux. Regions which lie above the line have more observed flux than is 
predicted, suggesting gyoresonance emission. Regions which lie below the line are non-physical and indicate 
more predicted optically thin bremsstrahlung emission than the total observed flux. Regions labelled in red 
have peak polarization brightness temperatures of Tb > 10 5 K, blue regions have 5 x 10 4 < Tb < 10 5 K, black 
regions have Tb < 5 x 10 4 K, and green regions lie above the solar limb. 


6. Discussion 
6.1. Full disk 

A dramatic feature of Figure[3]is that the EUV- 
predicted radio fluxes of the regions at the so¬ 
lar limb are all well in excess of their counter¬ 
parts in the radio image. This effect can also 
be seen in the left panel of Figure [5] where the 
limb regions are plotted as green points which 
all lie in the non-physical regime. The total 
predicted bremsstrahlung flux from the limb re¬ 
gions was 8.2 ± 0.3 sfu whereas only 2.3 ± 0.1 
sfu was observed. Correcting for this 5.9 ± 0.3 
sfu difference in the limb flux suggests that there 
should be 71.8 ± 0.3 sfu of optically thin coronal 
bremsstrahlung in the F 10.7 signal. Comparison to 
the observed 76.3±2.4 sfu variable component now 
suggests the presence of a small amount of gyrores- 
onance emission. Note that this is a conservative 
correction for the discrepancy at the limb because 
it does not account for any of the limb emission 


outside the main active regions. 

The chromosphere provides the optically thick 
background for both the radio and EUV obser¬ 
vations and therefore sets the height of the vis¬ 
ible solar limb, however this height is frequency 
dependent. The effective solar limb at 2.8 GHz 
is aro und 30" above the solar photosphere ( Gary : 
1996h , while the height of the effective solar limb at 


EUV wavelength s is only a few arcseconds above 
the photosphere ( Auchere et al.Ml998L Zhang ct al 


It is believed that the extra height of the ra¬ 
dio limb is due to cool filamentary chromospheric 
material (such as spicules) that extends into the 
solar corona and can provide extra opacity at the 
limb. We interpret the depressed radio signal from 
the limb regions as the occultation of emission 
originating behind the chromospheric limb. This 
has interesting implications for F 10.7 as an EUV 
proxy because the difference in limb altitude may 
cause the obscuration of a significant fraction of 
the solar emission in the radio which is visible in 
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the EUV. In this case, plasma that produced at 
least 7.7 ± 0.5% of the variable F 10.7 component 
and which was visible in EUV was not observed 
in the radio. Not only does this complicate com¬ 
parison of F 10.7 and EUV fluxes, it also causes 
an offset in time-series comparisons because EUV 
sources will become visible before corresponding 
F 10.7 sources rotate into view, and will remain vis¬ 
ible after F 10.7 sources rotate behind the limb. 


6.2. Individual regions 


Bremsstrahlung emission is usually weakly 
polarized and gyroresonance sources are often 
strongly polarized, therefore those regions which 
had highly polarized emission were suspected gy¬ 
roresonance sources. In Figure [5] the observed re¬ 
gions were categorized based on the peak polariza¬ 
tion brightness temperature. Regions with a peak 
polarization brightness temperature Tb > 10 5 K 
at 2.783 GHz were deemed to have strong polar¬ 
ization and all lie well above the bremsstrahlung 
line, therefore confirming that they have signif¬ 
icant gyroresonance contribution. Regions with 
5 x 10 4 < Tb < 10 5 K were classified as weakly po¬ 
larized, possibly containing gyroresonance sources, 
while regions with Tb<5x10 4 K were deemed 
to be insufficiently polarized and were unlikely to 
contain significant gyroresonance sources. These 
weakly polarized and unpolarized regions all have 
fluxes very close to the bremsstrahlung predic¬ 
tions, with only two (regions 5 and 9) showing 
significant observed radio excess. The total radio 
excess which is interpreted as the gyroresonance 
contribution (although this could also be due in 
part to the bremsstrahlung becoming optically 
thick) was 5.6 ± 0.2 sfu in the regions with strong 
polarization, 0.3 ± 0.2 sfu in the weakly polarized 
regions, and 0.3 ± 0.2 sfu in the unpolarized re¬ 
gions. More than 60% of the total gyroresonance 
emission originates in region 3, the largest disk 
active region. 


This analysis suggests that 6.2 ± 0.3 sfu or 
8.1 ± 0.5% of the variable F 10.7 signal recorded 
on 2011 December 9 resulted from gyroresonance 
emission. While this is a small percentage it is well 
above the precision of the Penticton F 10.7 mea¬ 
surements and could be sufficient to account for 
the known ~ 10 % density errors resulting from 
Fm 7 driven thermosphere models ([Bowman et al 


2008!) . It is important to note that this gyroreso¬ 


nance emission, if it were constant (or even a con¬ 
stant fraction of the variable component), would 
have little impact on the use of F 10.7 as an EUV 
proxy. However, since the gyroresonance contri¬ 
bution is dominated by the largest active region, 
we speculate that it is likely to be a much larger 
fraction of F 10.7 at times of high solar activity and 
this could affect the use of F 10.7 as an EUV proxy. 


6.3. Coronal Iron Abundance 


A straightforward result of this analysis is con¬ 
firmation that the iron abundance in the corona is 
iVpe/Vj-j = 1.26 x 10 -4 , an enhancement of about 
a factor of 4 over the photospheric value. Two 
independent results which depend on the coronal 
value of N-pe/N-B confirm this conclusion: the to¬ 
tal bremsstrahlung flux derived from the AIA data 
matches the coronal contribution to F 10.7 derived 
by subtracting the solar minimum flux from the 
measured F10.7 flux during the observation; and 
the fluxes of the 17 GHz active regions, which 
are known to be well calibrated by matching the 
NoRH full-disk flux to the calibrated NoRP patrol 
measurement, are consistent with the predictions 
from AIA data. The inferred DEM of H, to which 
the predicted bremsstrahlung radio flux is pro¬ 
portional, depends inversely on the assumed iron 
abundance: if Np e /Nn were to be photospheric, 
the F10.7 prediction based on the AIA data for 
this day would be 67.2 + (4 * 77.7) = 378 sfu, 
rather than the measure d 143.5 ± 1.2 sfu. T hese 
results are consistent witli lWhite et al.1 ( 2000 1 who 
found an iron abundance of Np e /Nn = 1.56 x 
10” 4 with approximate 20% errors by compar¬ 
ing radio observations with bremsstrahlung pre¬ 
dictions from EUV DEMs in a single active region. 
Note that our assumption that the solar minimum 
flux has no significant c oronal contr i bution dis- 
agrees with the model of Zirin et ah (1991), but 
Landi &; Chiuderi Dragol ( 20031 2008h carried out 
a careful comparison of the solar minimum radio 
spectrum with a DEM derived from UV and EUV 
data and in their results the F 10.7 solar minimum 
flux has only a small coronal contribution. 

Features other than coronal emission from ac¬ 
tive regions, including flares (which may be domi¬ 
nated by evaporated chromospheric material) and 
energetic particles, have shown differe nt abun¬ 
dances, ranging from 1.2 ( Meveil Il985l ) to 13.1 
( Reamesl 19991) times the photospheric level. Our 
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results suggests that an iron enhancement of 4 is 
generally appropriate for coronal active regions: 
we will pursue this result further in a future study. 


ence in the gyroresonance flux, still restricts the 
total gyroresonance to less than 15% of the vari¬ 
able F 10.7 component. 


6.4. Uncertainties 


7. Conclusion 


It should again be noted that the systematic 
uncertainties involved in this analysis are much 
greater than the statistical errors quoted above. 
Systematics such as the iron abundance and im¬ 
proper calibration of the VLA 20 dB attenuators 
cause constant offsets and therefore affect the over¬ 
all agreement of the disk fluxes (although this does 
not apply to the well-calibrated 17 GHz data). 
The effects of these constant offsets is minimized 
by normalizing each region to the local disk back¬ 
ground as described in section [5] However, based 
on the errors in coronal abundance studies as well 
as differences between the flux observed with the 
VLA and the official record measured at Pentic¬ 
ton, we expect both of these error sources to be 
on the order of 20%. 

Uncertainties in the DEM solutions could have 
spatially variable errors depending on the under¬ 
lying plasma parameters. We expect these er¬ 
rors to be on the order of 10% or less but it is 
difficult to quantify the extent to which devia¬ 
tions between the derived DEM and the ground 
truth plasma parameters change the results due 
to the non-linear influence of the DEM on the 
bremsstrahlung prediction. However, the compar¬ 
ison of the bremsstrahlung prediction to the 17 
GHz Nobeyama observation in Figure [5] shows no 
clear correlation of deviation with active region 
size. This suggests that there is no systematic 
bias and that any pixel scale statistical errors in 
the DEM are washed out when integrating over an 
entire active region. 

No attempt was made to account for the un¬ 
derestimation of gyroresonance emission because 
of the predicted bremsstrahlung emission origi¬ 
nating from below the optically thick gyroreso¬ 
nance layer. If the most extreme case is as¬ 
sumed, that all of the observed radio emission 
from the strongly polarized regions resulted from 
gyroresonance emission (i.e., all of the predicted 
bremsstrahlung emission occurred below the op¬ 
tically thick gyroresonance layer), then the total 
gyroresonance flux from these regions would be 
10.9±0.1 sfu. This is a generous upper limit which, 
while it does allow a possible factor of two differ- 


Understanding the sources of the solar F 10.7 
flux is important if it is to be used reliably as an 
EUV proxy in thermosphere/ionosphere models. 
By comparing a full disk VLA observation with 
the F 10.7 bremsstrahlung emission predicted from 
DEMs calculated with AIA images, we find that 
8.1 ± 0.5% of the variable F 10.7 flux on 2011 De¬ 
cember 9 can be attributed to the gyroresonance 
mechanism. This gyroresonance contribution does 
not directly correlate with solar EUV flux and 
therefore should be removed from F 10.7 when it is 
used as an EUV proxy. While this is a small frac¬ 
tion of the F 10.7 signal, it is commensurate with 
the density errors in current ionospheric modelling 
efforts. 

We also identify unexpected occultation due to 
the optically thick chromosphere of F 10.7 flux orig¬ 
inating from behind the solar limb. It appears that 
this effect could cause systematic errors in F 10.7 
time series at the active region scale (on the order 
of 10 sfu) at one day time scales. To our knowl¬ 
edge, this effect has not previously been consid¬ 
ered as a possible complication when comparing 
F 10.7 to EUV emission. For our observation, the 
on-disk gyroresonance excess ( 6 .2 ±0.3 sfu) almost 
perfectly cancels the off disk paucity (minimum of 
5.9 ± 0.3 sfu). There is no reason for these two 
effects to be correlated except in the general sense 
that both are likely to vary with the general level 
of activity, and consequently it seems unlikely that 
they will generally offset each other as well as they 
do here. 


The details of these results depend strongly on 
the coronal iron abundance which is inherent in 
the calculation of the DEM. By comparing the 
bremsstrahlung prediction with the coronal con¬ 
tribution to the F 10.7 measurement and with the 
Nobeyama 17 GHz images, we confirm a coro¬ 
nal iron abundance of N-p e /Nn = 1.26 x 10 -4 
(the standar d coronal abunda n ce in the CH I ANTI 


database; Dere et al. ] Il997t lLandi et al.1 l2013h 


which is used in the calculation of AIA temper¬ 
ature response functions. These results are sub¬ 
ject to various potential systematic error sources 
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which are difficult to quantify, but have estimated 
accuracies on the order of 20 %. 

No general statements about the effects of gy- 
roresonance and limb corrections on the F 10.7 in¬ 
dex can be made at this time because a single 
measurement is insufficient to characterize their 
temporal variability. Additionally, based on the 
untested systematic biases from the VLA calibra¬ 
tion complications and the confusion regarding the 
altitude of optically thick gyroresonance layers, 
caution should be taken when considering these 
results. The temporal variability of the gyrores¬ 
onance fraction and the effects of gyroresonance 
emission altitude will be addressed further in fu¬ 
ture studies for which the data have already been 
collected. 

Acknowledgements: S.W. would like to ac¬ 
knowledge valuable conversations with Ken Tap¬ 
ping. The National Radio Astronomy Observatory 
is a facility of the National Science Foundation 
operated under cooperative agreement by Associ¬ 
ated Universities, Inc. We thank the Nobeyama 
Solar Radio Observatory for providing ready ac¬ 
cess to their data. Data supplied courtesy of 
the SDO/HMI and SDO/AIA consortia. SDO 
is the first mission to be launched for NASA’s 
Living With a Star (LWS) Program. CHIANTI 
is a collaborative project involving George Ma¬ 
son University, the University of Michigan (USA) 
and the University of Cambridge (UK). This re¬ 
search has been made possible with funding from 
AFOSR FA9550-15-1-0014, NSF Career Award 
#1255024, New Mexico Space Grant (a NASA 
EPSCoR branch), and PAARE NSF:0849986. 

REFERENCES 

Asplund, M., Grevesse, N., Sauval, a. J., & Scott, 
P. 2009, Annual Review of Astronomy and As¬ 
trophysics, 47, 481 

Auchere, F., Boulade, S., Koutchnry, S., et al. 
1998, Astronomy & Astrophysics, 336L, 57 

Bastian, T. S., & Dulk, G. A. 1988, in Solar and 
Stellar Coronal Structure and Dynamics, ed. 
R. C. Altrock, 386-391 

Bhatnagar, V. P., & Mitra, A. P. 1966, Journal of 
the Atmospheric Sciences, 23, 233 


Boerner, P., Edwards, C., Lernen, J., et al. 2012, 
Sol. Phys., 275, 41 

Bouwer, S. D. 1992, Solar Physics, 142, 365 

Bowman, B. R., Kent Tobiska, W., Marcos, F. a., 
& Valladares, C. 2008, Journal of Atmospheric 
and Solar-Terrestrial Physics, 70, 774 

Chen, Y., Liu, L., & Wan, W. 2011, Journal of 
Geophysical Research, 116, A04304 

Covington, A. E. 1947, Nature, 159, 405 

—. 1948, Proceedings of the IRE, 36, 454 

—. 1951, Journal of the Royal Astronomical Soci¬ 
ety of Canada, 45, 15 

—. 1969, Journal of the Royal Astronomical Soci¬ 
ety of Canada, 63, 125 

Covington, A. E., Medd, W. J., Harvey, G. A., 
& Broten, N. W. 1955, Journal of the Royal 
Astronomical Society of Canada, 49, 235 

Craig, I. J. D., & Brown, J. C. 1976, Astronomy 
& Astrophysics, 49, 239 

Deng, L. H., Li, B., Zheng, Y. F., & Cheng, X. M. 

2013, New Astronomy, 23-24, 1 

Dere, K. P., Landi, E., Mason, H. E., Monsignori 
Fossi, B. C., & Young, P. R. 1997, Astronomy 
and Astrophysics Supplement Series, 125, 149 

Dere, K. P., Landi, E., Young, P. R., et al. 2009, 
Astronomy & Astrophysics, 498, 915 

Dudok de Wit, T., Bruinsma, S., & Shibasaki, K. 

2014, Journal of Space Weather and Space Cli¬ 
mate, 4, A06 

Dudok de Wit, T., Kretzschmar, M., Lilensten, 
J., & Woods, T. 2009, Geophysical Research 
Letters, 36, L10107 

Dulk, G. A. 1985, Annual Review, Astronomy 
&Astrophysics, 23, 169 

Felli, M., Lang, K. R., & Willson, R. F. 1981, The 
Astrophysical Journal, 247, 325 

Foukal, P. 1998, Geophysical Research Letters, 25, 
2909 


14 



Frohlich, C. 2009, Astronomy & Astrophysics, 
501, L27 

Furst, E., Hirth, W., & Lantos, P. 1979, Solar 
Physics, 63, 257 

Gary, D. E. 1996, in Astronomical Society of the 
Pacific Conference Series, ed. A. R. Taylor & 
J. M. Paredes, Vol. 93, San Francisco, CA, 387 
396 

Golub, L., & Pasaclioff, J. M. 2010, The Solar 
Corona, 2nd edn. (Cambridge University Press) 

Gopalswamy, N., White, S. M., & Kundu, M. R. 
1991, The Astrophysical Journal, 379, 366 

Greenstein, J. L., & Minkowski, R. 1953, The As¬ 
trophysical Journal, 118, 1 

Hannah, I. G., & Kontar, E. P. 2012, Astronomy 
& Astrophysics, 539, 146 

Henney, C. J., Toussaint, W. a., White, S. M., & 
Arge, C. N. 2012, Space Weather: The Inter¬ 
national Journal of Research and Applications, 
10, S02011 

Jacchia, L. G. 1971, Revised Static Models of the 
Thermosphere and Exosphere with Empirical 
Temperature Profiles, Tech, rep., Smithsonian 
Astrophysical Observatory 

Johnson, R. W. 2011, Astrophysics and Space Sci¬ 
ence, 332, 73 

Kundu, M. R. 1959, International Astronomical 
Union, 9, 222 

. 1965, Solar Radio Astronomy (Interscience 
Publishers) 

Landi, E., & Chiuderi Drago, F. 2003, The Astro- 
physical Journal, 589, 1054 

-. 2008, The Astrophysical Journal, 675, 1629 

Landi, E., Young, P. R., Dere, K. P., Del Zanna, 
G., & Mason, H. E. 2013, The Astrophysical 
Journal, 763, 86 

Lean, J. L., & Brueckner, G. E. 1989, The Astro- 
physical Journal, 337, 568 

Lenren, J. R., Boerner, P. F., Edwards, C. G., 
et al. 2012, Sol. Phys., 275, 17 


Martyn, D. F. 1948, Proceedings of the Royal So¬ 
ciety Series A Mathematical and Physical Sci¬ 
ences, 193, 44 

Maruyama, T. 2010, Journal of Geophysical Re¬ 
search, 115, A04306 

—. 2011, Journal of Geophysical Research, 116, 
A08303 

Meyer, J.-P. 1985, The Astrophysical Journal Sup¬ 
plement Series, 57, 173 

Nakajima, H., Sekiguchi, H., Sawa, M., et al. 
1985, Publications of the Astronomical Society 
of Japan, 37,163 

Nakajima, H., Nishio, M., Enonre, S., et al. 1994, 
IEEE Proceedings, 82, 705 

Parker, D. G., Ulrich, R. K., & Pap, J. M. 1998, 
Solar Physics, 177, 229 

Pesnell, W. D., Thompson, B. J., & Chamberlin, 
P. C. 2011, Solar Physics, 275, 3 

Piddington, J. H., & Minnett, H. C. 1951, The 
Commonwealth Scientific and Industrial Re¬ 
search Organization Australia, 4, 131 

Plowman, J., Kankelborg, C., & Martens, P. 2013, 
The Astrophysical Journal, 771, 2 

Poduval, B., DeForest, C. E., Schnrelz, J. T., & 
Pathak, S. 2013, The Astrophysical Journal, 
765, 144 

Reames, D. V. 1999, The Astrophysical Journal, 
518, 473 

Saint-Hilaire, P., Hurford, G. J., Keating, G., 
Bower, G. C., & Gutierrez-Kraybill, C. 2012, 
Solar Physics, 277, 431 

Scherrer, P. H., Schou, J., Bush, R. I., et al. 2011, 
Solar Physics, 275, 207 

Schmahl, E. J., & Kundu, M. R. 1995, Journal of 
Geophysical Research, 100, 19851 

Schmahl, E. J., & Kundu, M. R. 1998, in As¬ 
tronomy Society of the Pacific Conference Se¬ 
ries, ed. K. S. Balasubranraniam, J. Harvey, & 
D. Rabin, Vol. 140, San Francisco, CA, 387-399 

Schou, J., Scherrer, P. H., Bush, R. I., et al. 2011, 
Solar Physics, 275, 229 


15 



Sclhorst, C. L., Costa, J. E. R., Gimenez de Cas¬ 
tro, C. G., et al. 2014, The Astrophysical Jour¬ 
nal, 790, 134 

Svalgaard, L., & Hudson, H. S. 2010, Astronomy 
Society of the Pacific Conference Series, 428, 
325 

Swarup, G., Kakinuma, T., Covington, A. E., 
et al. 1963, The Astrophysical Journal, 137, 
1251 

Tapping, K. F. 1987, Journal of Geophysical Re¬ 
search, 92, 829 

—. 2013, Space Weather, 11, 394 

Tapping, K. F., Boteler, D., Charbonneau, P., 
et al. 2007, Solar Physics, 246, 309 

Tapping, K. F., Cameron, H. T., & Willis, A. G. 
2003, Solar Physics, 215, 357 

Tapping, K. F., & Clrarrois, D. P. 1994, Solar 
Physics, 150, 305 

Tapping, K. F., & DeTracey, B. 1990, Solar 
Physics, 127, 321 

Tapping, K. F., & Valdes, J. J. 2011, Solar 
Physics, 272, 337 

Tobiska, W. K., Bouwer, S. D., & Bowman, 
B. R. 2008, Journal of Atmospheric and Solar- 
Terrestrial Physics, 70, 803 

Vats, H. O., Deshpande, M. R., Shah, C. R., & 
Mehta, M. 1998, Solar Physics, 181, 351 

White, S. M. 1999, Solar Physics, 190, 309 

White, S. M., & Kundu, M. R. 1997, Solar Physics, 
174, 31 

White, S. M., Kundu, M. R., & Gopalswamy, N. 
1992, The Astrophysical Journal, 78, 599 

White, S. M., Thomas, R. J., Brosius, J., & 
Kundu, M. R. 2000, The Astrophysical Jour¬ 
nal, 534, 203 

Wild, J. P., Smerd, S. F., & Weiss, A. A. 1963, 
Annual Review, Astronomy & Astrophysics, 1, 
291 

Wilson, R. M., Rabin, D., & Moore, R. L. 1987, 
Solar Physics, 111, 279 


Zanna, G. D., Dwyer, B. O., & Mason, H. E. 2011, 
Astronomy & Astrophysics, 535, 1 

Zhang, J., Kundu, M. R., White, S. M., Dere, 
K. P., & Newmark, J. S. 2001, The Astrophys¬ 
ical Journal, 561, 396 

Zhang, J., White, S. M., & Kundu, M. R. 1998, 
The Astrophysical Journal, 504, L127 

Zirin, H., Baumert, B. M., & Hurford, G. J. 1991, 
The Astrophysical Journal, 370, 779 


This 2-column preprint was prepared with the AAS IATgX 
macros v5.2. 


16 



